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This paper investigates origin of the extra stability associated with particular values (magic 
numbers) of the total angular momentum of electrons in a quantum dot under strong magnetic 
field. The ground-state energy, distribution functions of density and angular momentum, and 
pair correlation function are calculated in the strong field limit by numerical diagonalization 
of the system containing up to seven electrons. It is shown that the composite fermion picture 
explains the small magic numbers well, while a simple geometrical picture does better as the 
magic number increases. Combination of these two pictures leads to identification of all the 
^sO • magic numbers. Relation of the magic-number states to the Wigner crystal and the fractional 

. quantum Hall state is discussed. 
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O ■ §1. Introduction 



The two-dimensional electron system with a small number N of electrons is realized as the 
quantum dot formed at a semiconductor interface.! Correspondingly there is a growing theoretical 
interest in the system. It has been found!!!!!!!!© that if a system is placed in a strong 
magnetic field, extra stability arises for a special set of angular momentum M. These values of 
M are called magic numbers. According to numerical diagonalization, the magic number states 
occur with an interval AM = N for systems with N < 5, and that these states have a polygonal 
pattern with iV apexes in the configuration of electrons.!! The magic number states prevail the 
ground-state phase diagram in the plane of magnetic field vs the confining potential for N = 5 
and 6^ On the other hand, it has been shown!! that many of the magic number states can be 
explained in terms of the composite fermion (CF) picture© for the fractional quantum Hall (FQH) 
states. The overlap is close to unity!© between a trial wave function based on the CF theory and 
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the exact ground state wave function in the presence of the Coulomb interaction. However, some 
magic number states do not fit into the interpretation of the CF theory. 

The purpose of this paper is to clarify the origin of the magic numbers on the basis of exactly 
derived results for the energy, and one- and two-body distribution functions. In addition to the 
previously known series with AM = N, we find another series with AM = N — 1 of the magic 
number states for N > 6. We provide simple geometrical interpretation on the origin of new series 
of magic numbers. Combining with previous interpretations, we come to unified understanding of 
all the magic numbers. 

The paper is organized as follows: In §2 we define the model which involves truncation to the 
lowest Landau level. Section 3 presents numerical results on the ground state energy as a function of 
the total angular momentum. The magic numbers are identified for systems with 6 and 7 electrons. 
The results on one- and two-body distribution functions are shown in §4 where the characteristic 
features of the magic-number states are made explicit. We provide in §5 a geometrical interpretation 
of the magic number states, and discuss the relation to the CF picture.^ In §6 we combine the 
geometrical and CF interpretations, and point out the range of validity for each interpretation. 
This reveals significance of magic-number states as showing crossover from an incipient fractional 
quantum Hall state to an incipient Wigner crystal. 

§2. Model and Choice of Basis 

We consider a model for a quantum dot in two dimensions with the Hamiltonian 

H = H + V 1 + V 2 . (2.1) 

Here Hq, V%, and V 2 are given by 

1 N h 

^o = ^ 7 E(-V, + -A i ) 2 , (2.2) 
2m* ~ 1 c 

i=i 

fi = imV5>| a , V 2 = -Y.r^- V (2-3) 

with ijj being the frequency of the harmonic oscillator, and z = x + iy. Other notations are the 
standard ones. The vector potential is chosen to be Aj = (jJiB/2,—X{B/2). We assume that the 
magnetic field B is so strong that all electrons are in the lowest Landau level (LLL), with spins 
being completely polarized. One-electron state is specified by its orbital angular momentum only. 
We ignore the spin degrees of freedom and mixing between the LLL and higher Landau levels. 
Hence the Zeeman term has been dropped in our model. 
The normalized single electron wave function in the LLL is given by 
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where m is the angular momentum and £ the magnetic length defined by £ = (fic/eB)^. The 
Hamiltonian in the second quantization is written as 



H = ^2 < m\Ho\m > a^dm + < m|Vi|m > a^ m a m 

+ X! < mi'm2\V2\m 3 m4 > a^a^a^a^, (2.5) 

mi 

where a m is the annihilation operator of an electron with m. The matrix elements are given by 

< m\H \m >= ^huj c , (2.6) 

< m\Vi\m > = Kn*uj 2 £ 2 (m + 1), (2.7) 
and the Coulomb matrix element is computed numerically from 

, e 2 . 2- d T{j-k+l)T(k + d+l) 
< mi m 2 \V 2 \m 3 m 4 >= (-) _ _ + ^ _ = _ = 

POO r*3> J>2 

x / dxx^F^j -k + l;d+l;-—) 1 F 1 (k + d+l;d+l;- — ). 
Jo 2 2 

Here \F\ is Kummer's hypergeometric function, and m\ = j — k,rri2 = k^m^ = k + d^m^ = j — k — d 
with d > 0. We discard in the following the constant kinetic term < m\Ho\m >= ^hu> c . The 
confinement parameter 7 is defined by 

m*uj 2 e 2 , . 

The length is scaled by £, and energy is scaled by e 2 /{eq£) from now on. 

Since the confining potential and the Coulomb interaction have the rotational symmetry, the 
total angular momentum M can be used to label the many body eigenstates. We introduce the 
distribution function of angular momentum by 

n m =< $ |<4om|$o >, (2.9) 

where |<£o > is the normalized ground state. Then the total angular momentum is given by 

00 

M = mn m- (2-10) 

m=0 

Next we introduce a field operator ip(z) truncated within the LLL: 

il>{z) =^> m (z)a m . (2.11) 

m 

Then the charge density n(z) is given by 

n{z) =< $ \il>i(z)il>{z)\$o >, (2.12) 



and the total number N by its integral over the whole space. Finally the pair (two-body) correlation 
function is defined as 

n(z, z') =< <f> \^(z)^(z')iJj(z'W(z)\$ > . (2.13) 
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§3. Two Series of Magic Numbers 

By using the Lanczos method to diagonalize the Hamiltonian numerically, we calculate the lowest 
energy and corresponding eigenvectors for each M and N. For given N and M, the one body angular 
momentum takes a value from to M — N(N — 1) /2. The upper limit less than M comes from the 
Pauli exclusion principle. In constructing the Hamiltonian matrix, we take all the necessary basis 
for m without truncation. We note the following relation: 

Vi\$ >=j(M + N)\$ >, (3.1) 

which follows from eqs.(2.7). Thus the wave function |<I>o > is actually independent of 7, and the 
numerical result for 7 = is sufficient to derive the energy for other cases of 7. 
Figure 1(a) shows the results for the ground-state energy vs M in the case of N = 6. A number 



Fig. 1. The ground-state energy vs the total angular momentum in the case of (a) N = 6 and (b) N = 7. The 
parameter 7 represents the strength of the harmonic confinement as explained in the text. 



of downward cusps appear clearly. Upon closer inspection one finds that a part of the cusps with 
M = 15, 21, 27, 33, 39, 45 are represented by the formula 

M = -N(N - 1) + Nk, (Jfe = 0,l,---) (3.2) 

with iV = 6. This set of states with the interval AM = N are called the series Sn in this 
paper. The other set of cusps appear for states with M = 25, 30, 35, 40, 45, 50 with another interval 
AM = N — 1 which is called the series Sn- The particular state with M = 45 belongs to both Sn 
and Sn. According to Fig. 1(a), M = 20 in the series Sn is not actually a magic number. With 
7 = 0, each magic-number state and the adjacent state with one more angular momentum often 
have the same energy. This degeneracy comes from the degrees of freedom associated with the 
center of mass motion as pointed out in ref]I|. 

The angular momentum of the absolute ground state is determined by competition between the 
Coulomb repulsion and the strength of the confinement . As 7 increases, the optimum angular 
momentum shifts to smaller values. In the case of 7 = 0.04 for example, the state with M = 45 
takes the absolute minimum of the energy. 
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In Fig. 1(b) we show the ground-state energy in the case of N = 7. We find clear downward cusps 
at the series S N . This includes M = 21, 33, 39, 45, 51, 57, 63 and 69. However, the number M = 27 
belonging to Sn is not a magic number. On the other hand, the cusps corresponding to the series 
Sn are not clear except for M = 28. This is in strong contrast with the case of N < 6. We note 
that a magic number state with M = 36 is seen in Fig. 1(b). This state belongs to neither the series 
Sn nor Sn • As we discuss later, the CF picture can interpret the magic number M = 36. 

§4. Distribution Functions 

4-1 Angular momentum 

In order to clarify the electronic property associated with the magic-number states, we calculate 
the distribution function n m of the angular momentum. In this section, we report mainly on results 
for N = 6. Figure 2 shows representative results. As M increases, two types of distribution appear 



Fig. 2. One-body angular momentum distribution function for (a) M = 33 and 35 and (b) M = 39 and 45. The 
distribution with double peaks as for M — 35 and 45 is called the type-D, while the one with a single peak as for 
M = 33 and 39 is called the type-S. 

by turns in the magic-number states. One is the distribution with double peaks at the origin as 
well as at a finite angular momentum. We call the double-peaked distribution the type-D hereafter. 
This is the case with M = 35 and 45. The other type is the distribution with a single peak which 
we call the type-S, as for M = 33 and 39. An important observation is that the magic number 
states with the type-D distribution are all in the series Sn, while those with the type-S distribution 
are all in the series Sn- These features are common to cases other than N = 6. 

It is instructive to interpret the result in terms of the Laughlin wave function . . . , zn) for 

finite N. It is given for general N by 



The odd integer p is related to the filling v = 1/p of the LLL in the limit of N — ► oo. In the case of 
finite N, the maximum one-body angular momentum m max is given by m max = pN(N — 1)/2. With 




(4.1) 
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p = 3 and N = 6, in particular, we obtain m max = 15 and the total angular momentum M is given 
by M = 3N(N — l)/2 = 45. We note that the numerically obtained results for M = 45 in Fig. 2(b) 
has the maximum 15 of the significantly occupied angular momentum in good correspondence to 
the Laughlin wave function with p = 3. 

4-2 Charge density 

The larger angular momentum corresponds to the wave function more extended from the origin. 
More explicitly we have the relation 

< m\r 2 \m >= 2£ 2 (m + 1). (4.2) 

We can therefore expect that n(z) of the type-D distribution should have double peaks both near 
the origin and near the edge of the dot, and that of the type-S distribution should have only a 
single peak near the edge. We have actually calculated the charge densities for these states with 
N = 6 and 7. Since the result n(z) depends only on r = \z\, the density is written as n(r) in the 
following. 

Figure 3 shows some exemplary results. It is found that the state with the type-S distribution 



Fig. 3. The charge density n(r) vs the distance r from the center of the dot in the case of M = 33, 35 and 45. 



(M = 33) indeed has a single peak near the edge of the quantum dot, and those with the type-D 
distribution (M = 35, 45) have double peaks. This feature is common to systems with different 
electron numbers. Since the Laughlin-type state with M = pN(N — l)/2 belongs simultaneously 
to series SV and SV, it is natural to expect that the charge distribution is also a superposition of 
the type-D and the type-S distributions. As a result, the distribution should have double peaks. 
The result in Fig. 3 shows that this is indeed the case with M = 45 and N = Q. 

It is seen from eq. ( |4.2D that the average density of electrons in the quantum dot decreases as 
the total angular mometum M increases. In §6 we use this fact for unified understanding of magic 
numbers. 
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4-3 Pair correlation function 

In order to see the correlated motion of electrons more closely we calculate the pair correlation 
function n(z, z') exactly. This quantity gives us a quantum analogue of a snapshot picture of 
electrons. More precisely, it gives the distribution of electrons on condition that one of the electrons 
is nailed down at z' . The pair correlation function for N = 5 or less has been calculated by 
Maksym,& who found a pattern corresponding to the polygon with N apexes for the series Sn- 
For N = 6 or more, however, we are not aware of exact results reported so far. Moreover nothing 
is known about the pair distribution for the series Sn- We define ro as the peak position of the 
charge density, which in the case of type-D distribution is to be taken at the position of the outer 
peak. Then we calculate n{z,z') numerically, setting \z'\ = ro. 

From Fig. 4 with N = 6, it is clear that electrons are arranged like a hexagon in the magic number 
state (M = 33) belonging to the series Sn, and like a pentagon in the magic number state (M = 35) 
of the series Sn- We recall that a Laughlin-type state with v = 1/p belongs simultaneously to the 



Fig. 4. Contour map of the pair correlation function n(z,z') for (a) M = 33 and (b) M — 35. The complex 
coordinate z' = roi corresponds to the center of the darkest region with ro = 2.199 in (a) and ro = 2.592 in (b). 
The interval of contours is 8.265 x 10~ 4 in (a) and is 6.472 x 1CT 4 in (b). 



series Sn and Sn- Then n(z, z') should be a superposition of pentagon- and hexagon-like patterns. 
To demonstrate this, the pair correlation function for p = 3 (M = 45) is calculated and is shown 
in Fig. 5 . It clearly confirms the expectation. 

Similar calculation is also performed for N = 7. Figure 6 shows the result for M = 39 which 
belongs to the series Sn- It shows the hexagon-type pattern as expected. On the contrary, we have 
confirmed that non-magic-number states such as M = 35 (N = 7) and M = 34 (N = 6) do not 
show a regular polygonal pattern. 

§5. Interpretation of Magic Numbers 

5.1 Geometrical interpretation 

As we have seen in the previous section, all the magic number states with N = 6 belong either to 
the series Sq or Sq. In this section, we discuss the relation between the total angular momentum and 
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Fig. 5. Contour map of the pair correlation function n(z, z') for M = 45 with z' — 2.915z. The interval of contours 
is 5.045 x 10" 4 . 



Fig. 6. Contour map of the pair correlation function n(z, z') for M = 39, N — 7 with z' = 2.442i. The interval of 
contours is 1.343 x 1CT 3 . 

the symmetry of the many-electron wave function. The argument extends the previous one which 
explains the series Sn- Let us consider the wave function $(zi, z%, • • • , z^) in the case where the 
coordinates correspond to apexes of a regular polygon: z n = Rexp(irufi) with n = 1, 2, ■ • • , N and 
(j) = 2tt/N. By applying exp(iM(f)) where M is the total angular momentum operator, we obtain 



The second equality follows from the antisymmetry of the fermion wave function. Then we get 



In order for the regular polygonal pattern with N apexes to be realized, $>(zi, zn) must be 

nonzero. This gives a selection rule on M as follows: 




Z N ) = $(zn,Z1,'" 



ZN-l) 



lexp(iM0) + (-l) N }<t>(z 1 ,z 2 ,--- 



zn) = 0. 



(5.1) 



M = N(j+ 1 -), (j = 0, !,-••) 



(5.2) 



if N is even and 



M = Nj, (i = 0, !,-••) 



(5.3) 



S 



if N is odd. The Pauli principle in addition requires M > N(N — l)/2. This selection rule leads 
to the series <Sjv- 

In order to examine the case of the polygonal pattern with N — 1 apexes for N electrons, we put 
zn = 0, z n = Rexp(in(fi') with n = 1, 2, • • • , N — 1 and $ = 2ir/(N — 1). Then we obtain 

^W^) + (-l)"" 1 ]*^!, 22, • ■ ■ , ^-1,0) = 0. (5.4) 
This leads to another selection rule: 

Af = (JV-l)j, (j = 0,l,---) (5.5) 

if N is even, and 

M = (iV-l)(j + i), (j = 0,l,---) (5.6) 

if N is odd. This selection rule leads to the series Sjy. 

These rules nicely explain the occurrence of the magic number states in the case of N < 6. 
However for larger N, the configuration with iV-apex pattern costs more energy than another apex 
pattern with one or more electrons in the interior. This is the reason why the magic number series 
57 is hardly seen in Fig. 1(b). The exceptional appearance of Sj is the magic number state M = 28 
with N = 7. The stability of the state, however, is better explained in terms of the CF picture as 
discussed below. 

Now we turn attention to the v = 1/3 Laughlin-type state which belongs to the series Sjy and 
5jv simultaneously. Since this state has superposition of the polygonal pattern with iV apexes and 
the one with N — 1 apexes, an extra stability is expected due to their resonant energy. 

5.2 Composite fermion picture 

In ref.p] an interpretation of magic numbers in terms of the CF picture is presented. In this 
picture the wave function is constructed in the following form: 

^cf(zi, ...,z n ) = Y[( Zi - z,-) 2 W (*i, • • • , zn), (5.7) 

i<j 

where ^o(zi, . . . , zn) is a wave function of free electrons, and q is a natural number. The Jastrow- 
type factor in eq.( ^.7[ ) is interpreted as binding a magnetic flux with strength 2q to each electron, 
hence the name of the CF. Although ^>q(zi, . . . , zn) is not restricted to the LLL, the projection 
operator V picks out only such component that belongs to the LLL. Thus the total angular mo- 
mentum M of ^cf(zi, ■ ■ ■ , z^r) is the sum of the part Mq associated with ^o^i, • • • , Zn) an d that 
coming from the Jastrow-type factor. Namely we obtain 

M = qN(N - 1) + M . (5.8) 

According to refj9| the magic number corresponds to such \Po(<zi, • • • , zn) that has compact occu- 
pation of each Landau level from the lowest possible angular momentum. In the case of iV < 5 the 
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CF picture explains all the magic numbers.© The CF magic numbers for N = 5 are the almost the 
same as the series S$ with additional ones M = 18 and 22 from the series £5. For N = 4 all magic 
numbers in S4 are also given by the CF picture, while the CF picture gives M = 12 as the only 
magic number in the series S4. The series Sn with the electron number N < 5 have been studied 
by other authorsBi'B as well. 

For the case of N = 6 and q = 1, the state with M = 45 has six electrons in the LLL. This 
state has Mq = 15. The other extreme case is the one where each of the six lowest Landau levels 
has only one electron, which we write as (1,1,1,1,1,1), resulting in Mq = — 15. If one excludes an 
occupation where a higher LL has more electrons than a lower one, the number of compact states 
is counted as 11 for N = 6. However, both occupations (3, 1, 1, 1) and (2, 2, 2) give M = 27, and 
both (4, 1, 1) and (3, 3) give M = 33. The set of the magic numbers with q = 1 are then 

M = 15, 21, 25, 27, 30, 33, 35, 39, 45. 

All these states belong to either Sq or Sq. With q = 2 the lowest magic number is M = 45 and the 
second lowest is 51. A nice feature of the CF picture is to reject M = 20 as a magic number in 
accordance with the numerical result in Fig. 1(a). On the other hand, the numerical result shows 
clear cusps at M = 40 and 50 in contradiction to the CF picture. 

Similar analysis is carried out for N = 7. The number of compact states is 15 for each q and the 
magic numbers with q = 1 are given by 

M = 21, 28, 33, 35, 36, 39, 41, 42, 43, 45, 48, 49, 51, 56, 63. 

By comparing with Fig. 2(b) we see that many of them correspond to cusps in the ground state 
energy. However, the numerical result shows clear cusps also at M = 57 and 69 both of which 
do not correspond to compact states. As N increases, the magic numbers predicted in the CF 
picture appear more densely than given in the series Sn and Sn- Some of them appear clearly in 
the numerical result, while others have almost no cusp. Unfortunately the relative stability among 
magic number states is not given in the CF picture. We remark that among the magic numbers 
which are not explained in the CF picture, M = 40 for N = 6 and M = 57 for N = 7 are both 
next to smaller magic numbers, 39 and 56. 

§6. Discussion and Conclusion 

We have seen that there are apparently conflicting interpretations of magic number states: the 
geometrical one which emphasizes the real space configuration, and the CF one which emphasizes 
the compact occupation of the one-body angular momentum. In the case of small M, the CF 
picture works better in rejecting M = 20 with N = 6, or M = 27 with N = 7. However, in the case 
of large M, the geometrical interpretation works better in explaining M = 40 and 50 with N = 6, 
and M = 57 and 69 with N = 7. In the intermediate values of M, many of the magic numbers are 
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common to both the geometrical and the CF pictures. By combining both pictures we can account 
for all the magic numbers. 

It should be noted that the smaller M corresponds to higher density of electrons in the macro- 
scopic limit. Thus the crossover in the effectiveness of the geometrical and the CF interpretations 
seems to reflect the transition from a FQH liquid at high density to the Wigner solid at low density 
in the macroscopic limit. The success of both interpretations in the intermediate range of M is 
intriguing. This suggests a smooth crossover for finite N from the incipient FQH liquid to the 
incipient Wigner solid with decreasing density. 

In this connection we remind the early wor k0) which argues that particular densities for the stable 
FQH states also stabilize the Wigner crystal with large zero-point motion and the ring exchange. 
It is thus not surprising that some of the electronic states in the quantum dot connect both to the 
FQH state and the Wigner crystal. Experimentally, a particular advantage of the quantum dot 
is the wide range of controllable electron numbers N. Thus one may observe how the electronic 
states change from the atomic type to the macroscopic type as ./V is increased. This feasibility is 
in strong contrast to genuine atoms. 

In summary, we have investigated the quantum dot system in a strong magnetic field by numerical 
diagonalization of the Hamiltonian. We have provided unified understanding of all the magic 
numbers in the total angular momentum. 
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